function Q = cal_Q(N,x_max)
h = x_max/N;
Q = zeros(N,N);
Q(1,N) = func_a(-0.5*h);
Q(1,1) = -func_a(-0.5*h)-func_a(0.5*h);
Q(1,2) = func_a(0.5*h);
for i = 2 : N - 1
    Q(i,i-1) = func_a((i-1-0.5)*h);
    Q(i,i) = -func_a((i-1-0.5)*h)-func_a((i-1+0.5)*h);
    Q(i,i+1) = func_a((i-1+0.5)*h);
end
Q(N,N-1) = func_a((N-1-0.5)*h);
Q(N,N) = -func_a((N-1-0.5)*h)-func_a((N-1+0.5)*h);
Q(N,1) = func_a((N-1+0.5)*h);
Q = Q/h^2;